arXiv:1507.07691vl [cond-mat.soft] 28Jul2015 


Current-voltage characteristics of seven-helix proteins from a cubic array of amino 

acids 

Eleonora Alfmitc0 

Dipartimento di Ingegneria dell’Innovazione, Universita del Salento, via Monteroni, 1-73100 Lecce, Italy 

Lino ReggianQ 

Dipartimento di Matematica e Fisica, "Ennio de Giorgi”, 

Universita del Salento, via Monteroni, 1-73100 Lecce, Italy 
(Dated: July 29, 2015) 

The electrical properties of a set of seven-helix transmembrane proteins, whose space arrangement 
(3d structure) is known, are investigated by using regular arrays of the amino acids. These structures, 
specifically cubes, have topological features similar to those shown by the chosen proteins. The 
theoretical results show a good agreement between the predicted current-voltage characteristics 
obtained from a cubic array and those obtained from a detailed 3D structure. The agreement is 
confirmed by available experiments on bacteriorhodopsin. Furthermore, all the analyzed proteins are 
found to share the same critical behaviour of the voltage-dependent conductance and of its variance. 

In particular, the cubic arrangement evidences a short plateau of the excess conductance and its 
variance at high voltages. The results of the present investigation show the possibility to predict the 
I-V characteristics of multiple-protein sample even in the absence of a detailed knowledge of their 
3D structure. 


I. INTRODUCTION 

The interest in the structure and function of a given 
protein is crossing the boundaries of biology, joining to 
electronics, applied physics and engineering. This rel¬ 
evant example of cross-fertilization is due to the quest 
of finding new and more efficient materials to produce 
green and renewable fuels 0 0, to reduce wasting [:’> . 
to design highly specific prostheses for medicine fH 0 , to 
insure human and animal safety 00. Proteins perform 
in vivo activities that, when introduced in electronic de¬ 
vices, should allow astonishing performances 0 , 0 . In 
particular, many transmembrane proteins act as highly 
specialized sensors, with levels of sensitivity and selec¬ 
tivity inconceivable for any artifact. The use of pro¬ 
teins in nano-bio-devices is, therefore, the topic of an 
emerging branch of electronics, recently introduced as 
proteotronics jioj . Proteotronics has the final objective 
of producing user-friendly, versatile devices for the daily 
life so as point-of-care testings and low invasive devices 
for medicine. As a matter of fact, the final device is con¬ 
ceived to be small and able to produce a clear response. 
This is allowed by the use of quite small active elements 
(protein samples), and an appropriate miniaturized elec¬ 
tronic interface able to convert a biochemical response 
into an electrical signal. Recent experiments have shown 
how the capture mechanism of a ligand (small molecules, 
odours, photons, etc.) has a solid conjecture concern¬ 
ing the expected protein electrical responses SEMI. 
The microscopic interpretation is actually based on the 
knowledge of the detailed tertiary structure of the single 
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protein, i.e. of the space location of its amino acids. How¬ 
ever, from one hand, such a knowledge is largely incom¬ 
plete and even in the few favorable cases is still a matter 
of active research far from a general consensus in the sci¬ 
entific community M- From another hand, this knowl¬ 
edge is an essential input for modeling the charge trans¬ 
port of protein-based devices ill 1 1 ldl 1 16H 19l ] . Indeed, the 
specific protein topology produces an electrical response 
very peculiar of the given protein IE EH BO. Thus, 
when this input is not known from a credited source, an 
approximate structure template able to reproduce the es¬ 
sential features of the electrical response, can be useful 
to model the protein electrical properties. In particu¬ 
lar, we conjecture that for proteins with similar size and 
shape (globular, fibrous, transmembrane, etc.) the gen¬ 
eral trend of the electrical response or, at least, of the 
current-voltage (I-V) characteristic, can be modeled by 
using a single template. 

The aim of this paper is to test the above conjecture. 
To this purpose, we have investigated three proteins be¬ 
longing to the transmembrane family whose model of the 
tertiary structure (hereafter simply indicated as the 3D 
structure) is known F° r one of them, the bac¬ 

teriorhodopsin (bR), the TV characteristic is also avail¬ 
able from experiments full [2~iT| . For these proteins the I-V 
characteristic has been calculated using an impedance 
network protein analogue (INPA) [Icij, [3 ■ The results 
obtained from the 3D structure have been compared with 
those obtained by using a regular structure, here taken as 
a cube, and proved to capture the main features of the 
protein electrical properties. The simple cube is built 
with a number of nodes as close as possible to the num¬ 
ber of amino acids given by the public data bank (PDB) 
3- At present only the response corresponding to the 
protein native state will be analyzed. The presence of 
the ligand modifies the protein structure and also its free 
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estimated value is $ = 216meV [T^j. 

The transition probability for tunneling assumes the 



FIG. 1. Cartoon (not to scale) of the cubeS and of an amino 
acid helix. Both structures have the same number of nodes. 

energy [l9[ and this is beyond the aims of this paper. 

The paper is organized as follows. In section [TT] the 
theoretical model is briefly recalled. In section IIIII the 
given proteins and the analogue cubic structures are in¬ 
troduced. Section EV| reports the main results of this 
study and final conclusions are summarized in Sec. [Vj 


where Vij is the potential drop between the couple of 

i,j amino acids, a = 2 Tj and m is the electron 
effective mass, here taken the same of the bare value. 
When the FN tunneling is applicable (eVij > <E>), the 
maximal resistivity decreases, following the rule: 

p( V) = PMAx{-rr—) + Pmin (1- 77 —) (3) 

where pmax = 4 x 10 13 flA. is the maximum value of 
resistivity taken to fit the TV characteristic at the lowest 
bias, and p m in = 4 x 10 5 flA the minimum resistivity 
value taken to fit the I-V characteristic at the highest 
voltages. 


II. MODEL 

The microscopic approach makes use of the INPA 
model that describes the single protein electrical prop¬ 
erties starting from the protein topology. The model was 
already detailed in previous papers [13, |22|, [23| In partic¬ 
ular, the amino acid space position, as given by the 3D 
protein structure, is used to build up a graph, where each 
node corresponds to an amino acid. Two nodes are then 
linked if their distance is below a cut-off value, Rc■ Fi¬ 
nally, the links are converted into elemental impedances, 
with the ensemble of impedances leading to an impedance 
network specific to describe the electrical properties of 
the single proteins. By appropriate scaling, the result 
of the single protein are compared with available experi¬ 
ments carried out on macroscopic samples. 

In the present case, since the objective is the static I-V 
characteristic, the impedances reduce to resistances. The 
observed current response in monolayers of bR validates 
a tunneling-like mechanism of charge transport, which 
has been successfully simulated by using a Monte Carlo 
technique [H, HH, [26|. The model describes the observed 
transition between a quasi-linear (Ohmic) regime at low 
bias, ruled by the direct tunneling mechanism (DT), and 
a super-Ohmic regime at high bias, mostly due to the in¬ 
jection or Fowler-Nordheim (FN) tunneling mechanism. 
The critical value of the bias controlling the cross-over 
between the regimes, Vic, depends on the network struc¬ 
ture and corresponds to a microscopic barrier height, $, 
taken to be independent from the specific structure. The 


III. MATERIALS 

We have analyzed three proteins, namely: bacteri- 
orhodopsin (bR), rat olfactory receptor (OR) 17, and 
bovine rhodopsin (BR), all in their native state, i.e. the 
state pertaining to the protein in the absence of a specific 
external stimulus (ground state) [lj|. As a consequence 
of an external stimulus, these proteins may change their 
native state into an active state. Bacteriorhodopsin is 
analyzed by using its 3D structure, as given by public 
PDB, in particular the 2NTU entry[l5j]. This structure 
consists of 222 amino acids, exhibits seven helices ar¬ 
ranged to form a loop and has a size of about 5,5 nm. 
Cube6 , the corresponding regular structure, consists of 6 
nodes for edge, making a total of 216 nodes. The dis¬ 
tance between contiguous nodes, the gen, is chosen of 
5,5 A which corresponds to a relative maximum of the 
bR amino acid distance distribution. With an analogous 
procedure the rat OR 17, whose 3D structure consists of 
327 amino acids [20j |. is compared with the cube7, the 
corresponding regular structure that consists of 7 nodes 
for edge, making a total of 343 nodes with gen= 5,5 A. In 
both these cases, the value of the interaction radius, Rc, 
has been chosen of 6 A as reported in previous papers 
[27l l2Sl . The same is done for BR, whose 3D structure 
consists of 348 amino acids, and the comparison is carried 
out with the cube 7. 

A cartoon of cubeS with 3 nodes for edge, and an a- 
helix with the same number of amino acids (27) is shown 
in Fig. [I] 
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FIG. 2. I-V characteristics of bacteriorhodopsin (empty cir¬ 
cles) and cube6 (continuous red line), Rc= 6A. The bold black 
line indicates the experimental data [U], color online. 



FIG. 3. I-V characteristics of rat OR 17 (empty squares) and 
cube 7 (continuous red line), Rc = 6A, color online. 


IV. RESULTS 

The calculated TV characteristics of cube6 and bR and 
of cubel and rat OR-I7 are reported in Fig.[2]and Fig. [3] 
respectively. 

In Fig. [2] available experimental data [3, [2lJ are also 
shown. To our knowledge, the measurement of the I-V 
characteristics on rat OR-I7 samples do not have been 
performed, yet. In both cases calculations by using the 
3D structure models for these proteins is well reproduced 
by the corresponding cubes, as well as the DT-FN tran¬ 
sition. The specific critical bias values are summarized 
in Tab. |H 

The normalized excess conductance, T = 
((g) — go) / gmax, is here taken as the order param¬ 
eter that describes the associated phase transition 
between the DT and FN regimes. With this definition, T 
assumes a value in the interval (0,1). The normalization 
is performed by using go, the lowest conductance value 


(taken at the bias value of 1CV 6 V ) and g m ax , the 
highest conductance value (taken at the bias value of 
10 4 V) which corresponds to the asymptotic Ohmic 
behaviour when all the resistivities of the network take 
the minimum value. 

Close to the phase transition, the order parameter cal¬ 
culated for bR and rat OR 17 exhibits a power-law be¬ 
haviour, when plotted as a function of the reduced bias, 
e = (V — Vc)/Vc- These results are reported in Fig. [4] 
As is usual in finite-size systems and due to the large fluc¬ 
tuations, rounding and shifting effects [29| do not allow 
to exactly pinpoint the phase transition. Furthermore, 
since the DT regime overlaps with the FN regime near 
the transition region, there can be identified two differ¬ 
ent ’’critical” bias values, Vc and Vq [26j]. The critical 
bias Vc signals the beginning of the FN regime, where 
the system starts to deviate from the Ohmic behaviour 
[26}, while Vq , the Ginzburg voltage [H, J30j, indicates 
the consolidation of the FN regime. 

A further signal of the existence of a universal 
behaviour is found by analyzing the scaled variance 
((g 2 ) — (s) 2 ) Lis) 2 - As just observed by analyzing differ¬ 
ent proteins |26|, this quantity shows a power law with 
a critical exponent close to 3. Results are reported in 
Fig. [5] From the figures above it also emerges that, better 
than the current, T reveals the main differences between 
the 3D and the cube structure. Indeed, while the former 
shows a continuous increase after the DT-FN transition, 
the latter exhibits a plateau at intermediate bias. We 
suggest that this plateau is associated with the better 
regularity of the cube with respect to the 3D structure. 
As a matter of fact, by choosing gen = 5, 5 A, and Rc= 
6 A, all the links have the same length and the pattern of 
potential barriers has itself a cubic shape. Conversely, in 
3D structures the branches have a wide range of lengths 
and resistances. Therefore, as the bias increases (over 
the FN regime), at least one link able to get the lowest 
resistivity is available. Thus the transition toward the 
final Ohmic regime takes place continuously. Conversely, 
in the chosen cubic structures, the same bias increase 
produces a stalemate: all the branches have the same re¬ 
sistance value and none of them is able to prevail and gain 
the minimal resistivity. As a consequence, a kind of crit¬ 
ical slowing down is observed at about 3 Vc- Finally, this 
condition is broken and all the branches achieve the min¬ 
imal resistivity, and a percolation-like transition is ob¬ 
served. This analysis is represented in figured where the 
scaled resistance, R/Rmax > is reported both for cube6 
and cube7. A similar result was obtained by simulating 
colossal magnetoresistance in a regular 2D spin lattice, 
with electrons jumping potential barriers placed in all the 
sites [3l| . Also in that case, the signature of a percolation 
transition was identified. 


As a third case of interest, the bovine rhodopsin (BR) has 
been analyzed. This protein 3D structure, entry 1U19 of 
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Cube6 

bR 

Cube7 

OR 17 

BR 

V c 

2,1 

1,7 

2,6 

3,8 

10 

links 

540 

692 

882 

1005 

1034 

nodes 

216 

222 

343 

327 

348 


TABLE I. The critical bias Vc, links and nodes for the ana¬ 
lyzed structures. Voltages are given in Volts, Rc= 6A. 



Cube7 

OR 17 

BR 

Vc 

2,6 

1,4 

2,6 

links 

882 

1339 

1417 

nodes 

327 

343 

348 


TABLE II. The critical bias Vc, links and nodes for the ana¬ 
lyzed structures. Voltages are given in Volts, Rc=7A. 


PDB, consists of 348 amino acids, similar to the rat OR 
17. Therefore, its I-V characteristics should be also de¬ 
scribed by the cube 7. This protein entry gives a structure 
a little bit more dilated than that of rat OR 17 and the 
I-V calculated by using Rc = 6 A exhibits the value 
Vc = 10 V, significantly larger when compared with the 
Vc = 2,6 V value obtained with cube7, Rc = 6 A (see 
Tab. [IJ. On the other hand, by increasing the value of Rc 
to 7 A, both the 3D structure and the cubel give the 
same value for Vc, similar to that calculated for rat OR 17 
by using Rc = 7 A. Notice that the results obtained with 
the cube structure are quite insensitive to the change of 
the Rc value inside the range (gen, \[2gen). Results are 
reported in Fig. |7] and Table El 

In summary, we have found that proteins with a similar 
number of amino acids exhibit similar I-V characteristics, 
well reproduced by the cubic structure with gen=5,5 A . 
The choice of a value of Rc = 6 or 7 A, is irrelevant for 



FIG. 4. a) Excess conductance for cube6 (full circles) and 
bR (open circles), with Rc= 6 A; b) Excess conductance for 
cube7 (full squares) and OR-I7 (open squares), with Rc= 6 

A 



m=3.35(0.09) m=3.37(0.04) 

FIG. 5. Scaled conductance variance, a) bR, empty circles, 
and cube6, full circles; b) rat OR 17, empty squares, and cube7, 
full squares. 



FIG. 6. Normalized resistances of cube6 and cube7. 


the cubic structure, which produces similar responses, 
but is noteworthy for the 3D structures. 


V. CONCLUSIONS 

The I-V characteristics of three seven-helix proteins, 
taken as a prototypes of this family, have been investi¬ 
gated within an INPA model. A comparative analysis 
between the 3D structure and a cubic arrangement of 
the amino acids is then carried out. The analysis shows 
a general agreement between the results of the 3D struc¬ 
tures and those of the cubic arrays that contain a similar 
number of amino acids. In particular, the strong non¬ 
linearity of the I-V characteristic is well reproduced, to- 
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FIG. 7. I-V characteristics of BR (circles), rat OR 17 (trian¬ 
gles) and cube! (continuous red line), Rc = 7A, color online. 


gether with some interesting peculiarities related to the 
phase transition between a direct and injection tunneling 
regime of charge transport. Accordingly, a cubic protein 
structure in the form presented here should provide real¬ 
istic expectations of the I-V characteristics ascribed to a 
given protein even in the absence of a detailed knowledge 
of its 3D structure, thus being of valuable help to predict 
macroscopic electrical properties. A further step towards 
the advancing of proteotronics. 
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